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Abstract 



00 ' This paper will examine a model with many agents, each of whom has a different 

belief about the dynamics of a risky asset. The agents are Bayesian and so learn 
about the asset over time. All agents are assumed to have a finite (but random) 
^ | ■ lifetime. When an agent dies, he passes his wealth (but not his knowledge) onto his 

heir. As a result, the agents never become sure of the dynamics of the risky asset. 
We derive expressions for the stock price and riskless rate. We then use numerical 
examples to exhibit their behaviour. 



1 Introduction 

This paper will look at a model of agents with heterogeneous beliefs. We assume 
that there is a single risky asset that produces a dividend process. Agents are 

On 

o 



unsure of the dynamics of the dividend process. Specifically, they do not know one 
of the parameters that governs its dynamics. Agents therefore form beliefs about 
this parameter and update these over time. To avoid agents eventually determining 
the true value of the parameter, we assume t hat agents are fin i te live d. 

The paper will build on previous work of Brown Sz Rogers] ( 20091 ). That paper 



explained the general theory of how to incorporate heterogeneous beliefs into a dy- 
namic equilibrium model. However, in the case in which the agents were Bayesian, 
it was seen that the agents would eventually determine the true drift of the divi- 
dend process. The purpose of this paper is therefore to investigate a model in which 
there is a non-trivial steady state. This is done through the assumption that the 
different agents are in fact dynasties. Each member of the dynasty has a finite but 
random lifetime and when that member dies, he will pass on his wealth, but not 
his knowledge, to his heir. The paper will explain how to construct and solve this 
model and will lead to a stationary distribution for the stock price. 
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As in iBrown fc Eogersl tod ), we assume that there is a single risky asset which 
pays a dividend continuously in time. In addition there is a riskless asset in zero net 
supply. The dividend process of the stock is now assumed to be a quadratic function 
of an Ornstein-Uhlenbeck (OU) process. All the agents know all the parameters of 
the OU process except the mean to which it reverts. All the agents observe the OU 
process as it evolves and so as time progresses they update their beliefs about the 
unknown parameter. However, since they are finite lived, they will never find its 
true value. 

The model described is quite si mple, yet already t here is enough to make the 
asset pricing non-trivial. Just as in lBrown &: Roger 1 tod ), the agents maximise 
their expected utilities subject to their budget constraints and we use these optimi- 
sation problems to derive a state price density. Using this state price density we can 
then price the risky asset as the net present value of future dividends. Comparative 
statics allow us to see how the stock price depends on the parameters of our model. 
We also produce a volatility surface for the stock, which behaves very reasonably. 

The structure of the paper is as follows. We give a brief literature review below. 
Section[2]introduces the model and solves the equilibrium to determine a state price 
density. Section [3] then uses this state price density to calculate the prices of the 
stock and bond; these calculations are non-trivial. Section H] looks at comparative 
statics of the model and section [5] concludes. 



1.1 Literature Review 



There is a large literature on h eterogeneous bel iefs, which h a s been discu s sed in 
detail i nlBrown fc Eogersl d200flh . Work i nclude s lKur3 (l2008bh iKurzl (|l994l\. iKurd 



(1 19971). iKurz & Mqtolesef 



2006). IKurzl t008a). 



Kurz et al 



(200 



E 



Fan 




Harrison fc Kreosi (11978T). iMorrisI Ill99ffl. Wu fc Quo! (I200.t). IWu fcGuc ( 
Harris k Ravivl (|l993h . lKandel Pearson! (jl995h . lBuraschi fc Jiltsovl (|2006h . |jouini Napp 
(|2007h . Closer to the work presented here are the papers that assume that there is 
a parameter of the economy that is unknown to the agents. We briefly review such 
mo dels he re. 

Basakl hood ) considers a two-agent model in which each agent receives an en- 
dowment process. There is also an extraneous process that agents believe may 
effect the economy. The endowment process and all its parameters are observed. 
The extraneous process is observed, but the parameters of the stochastic differen- 
tial equation (SDE) that drives it are not known to the agents. They form beliefs 
about the drift term in this SDE and update their beliefs in a Bayesian manner. 
The paper analyses this problem and derives quantities such as the consumption, 
the state price density and riskless rate. Basak also explains how to generalise the 
mo del to multip le agents and multiple extraneous processes. 



Basakl (|2005l ) also considers a model with two agents, who each receive an en- 
dowment process. The aggregate endowment process is observed by the agents. 
They also observe its volatility, but not its drift; they use filtering to determine 
this drift. There is assumed to be a bond and risky security, both in zero net 
supply. Again, agents do not know the drift of the stock price. Agents maximise 
the expected utility of consumption. He then solves for the equilibrium and uses 
it to derive interest rates and perceived market risk of the agents. He also gives 
a number of generalisations to the model. For example, he considers the case in 
which there is a process which does not directly affect the asset prices. However, 
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each agent thinks that this process does affect the dynamics of the asset prices and 
so this changes the equilibrium. He also looks at the case of multiple agents and 
again derives the riskless rate and perceived market prices of risk. The final part 
of his paper looks at further extensions to his model; for example, he explores a 
monetary model in which there is a money supply that is stochastic and agents 
dis agree on its drift. 

Gallmever Sz Hollifield ( 20081 ) have considered the effects of adding a short-sale 
constraint to a model with heterogeneous beliefs. They consider a model with 
two agents. These agents are unsure about the drift of the output process of the 
economy. They start with initial beliefs about the drift and use filtering to update 
these. The agent who is initially more pessimistic is assumed to have logarithmic 
utility and a short sale constraint. The optimistic agent is assumed to have general 
CRRA utility and does not have a short sale constraint. The authors examine 
this model and derive expressions for the state price densities, stock price and 
consumption. In particular, they examine the effects of the imposition of the short 
sale constraint o n the stock price . 

The paper of lZapaterol (| 19981 ) considers a model in which there is an aggregate 
endowment process that obeys an SDE driven by two independent Brownian mo- 
tions. The constant drift of the process is unknown to the agents. There are 2 
groups of agents and they each have a different Gaussian prior for this drift. Zap- 
atero also considers the case in which as well as observing the endowment process, 
the agents also see a signal, which again is driven by the two Brownian motions, 
but has unknown drift. Again, agents have prior beliefs about this drift, which they 
update. He derives an equilibrium and shows that volatility of the interest rate is 
higher in an economy with the additional information source. 

ILl ( 2007 ) considers a model with 2 groups of agents. There is a dividend process 
which obeys some SDE, but the drift of this SDE is unknown. The drift can satisfy 
one of two different SDEs. Each group of agents attaches a different probability 
to the drift obeying the two different SDEs. They update this probability as they 
observe more data. Agents are assumed to have log utility and Li derives the 
stock price, wealth and consumption of agents in this model. He also analyses the 
volatility of the stock price. 

Turning to the Bayesian learning side of our story, we remark that there is 
an extensive literature on Bayesian learning in finance and eco nomics in which 
agents update their beliefs as the y observe data. Work in cludes lHautsch fc Hess 
(120041 1 . iKandel fc Pearson! dX995h . ISchinkel et all tooi ). iKalai Lehreri (|l993h 
eac h of whom use s this Bayesian learning in quite different setups. For exam- 
ple, ISchinkel et al. (120021 ) apply Bayesian learning to n competitive firms who set 
prices but do not know the demand function. They observe demand at each step 
and use this to update their posterior belief for the state of the world, which then 
impacts their p e rceive d demand function. The authors show that prices converge. 
Kalai Lehrer ()l993h applies Bayesian learning to an n-person game in which 



agents do not know the payoff matrices of their competitors. Th ey show that the 



Hautsch fc Hess 



equilib rium will approach the Nash equilibrium of the system. 
( 20041 ) apply Bayesian learning to explain why more precise data has a larger im- 
pact on market prices. They test this by looking at the behaviour of T-bond futures 

when unemployment da ta is announced. 

Closer to our work, Guidolin Sz Timmermann ( 200ll ) look at a discrete time 
model in which the dividend process can have one of two different growth rates over 
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each time period and the probability of each growth rate is unknown to the agents. 
The agents are learning, so update their estimate for the unknown probability 
at each time step. In order to avoid the problem of agents discovering the true 
probability, they also consider agents who only look at a rolling window of data. 



2 The Model 



The setup of our model is similar to Brown Sz Rogers] ( 20091 ). There is a single 



productive asset, which we refer to as the stock, which pays dividends continuously 
in time. The dividend at time t is St. The dividend process is assumed to be a 
quadratic function of a stationary Ornstein Uhlenbeck (OU) process. 

Since we are interested in obtaining a stationary distribution for the stock 
pr ice, the construct i on of the probability space requires slightly more care than 
Brown Rogers (2009). Let $7 denote the sample space. We set f2 = 



m 

the space of continuous functions from R to R. Let X t (oj) = Lo(t) denote the 
canonical process. Furthermore, let Tt = cr(X s : — oo < s < t). 

As before, the reference measure is denoted by Po- We assume that under this 
measure X is a stationary OL0 process which reverts to mean zero and has reversion 
rate A. 

Next, we define: 

W t = X t -X + I \X s ds (2.1) 



o 



for all i € K. Since X is an OU process, we observe that the process (Wt)t>o is a 
standard Brownian motion^. 



2.1 The dividend process 

We now define the dividend process by: 

St = a + a\X t + a%Xf 

for some constants ao,ai,a2- 

The simplest non-trivial setup is that in which 00 = ^2 = 0, in which case the 
dividend process will simply be an OU process. However, choosing such values of 
ao and ci2 means that there is a positive probability that the dividend process will 
become negative, which is unrealistic. To overcome this problem, the constants can 
be chosen so that ao > af/4a2, in which case the dividend process will always be 
non- negative. Furthermore, it will transpire that considering the case in which the 
dividend process is a quadratic function of X is no more difficult than the case in 
which 5 is simply a scaling of 

*An Ornstein Uhlenbeck process which reverts to mean a' with reversion rate A satisfies the SDE 
dX t = dWt — A(a' — X t )dt where W is a standard Brownian motion under the reference measure. 

2 It will transpire that we are only interested in the increments of W; thus it does not matter that Wo 
is known before time 0. 

3 The case in which S is a quadratic function of X is slightly more complicated, since two different 
values of X can give the same value of 5. Hence, a(X s : to < s < t) ^ a(6 s : to < s < t). Thus, we must 
assume that the agents observe the process X , rather than just observing the process S. 
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2.2 The Agents 

In our model there are N agents at all times. We assume that each person has a 
random lifetime. When this person dies, their wealth is immediately passed onto 
their (ignorant) child. Thus we are viewing each agent as a dynasty rather than a 
persorQ. 

Formally, there exist times (Tl)k^z which are the jump times of a stationary 
renewal process. At each of these times Tfc, agent i will die and be replaced by 
his child. Thus, the wealth of the agent will be maintained, but their beliefs will 
not; the child will start with his own ignorant beliefs which will not depend on any 
historical data. 

Turning now to the beliefs of the agents, first recall that, under the reference 
measure, (Xt)t^R is an OU process with zero mean. However, under the true 
measure, X will revert to level a, which will not necessarily be zero. The agents do 
not know this level. They will use Bayesian updating to deduce it. 

We need to determine the measure that each agent works under. First note that 
if we restrict to the time interval [s,t], we may define a new measure by: 

^ = exp (\a(W t - W s ) - \(Xaf(t - s)\ (2.2) 

It follows from the Cameron-Martin-Girsanov theorem^ that a standard Brownian 
motion under Po becomes a Brownian motion with drift Xa under P a . Formally, 
W r = W r + Xar, for s < r < t where W is ci standard Browniaxi motion under IPq,. 
Thus, 

dX t = dW t + X(a - X t )dt 

so we see that, under P , X is an OU process which reverts to mean a. 

Since agents do not know a, the beliefs of each agent simply consist of their 
distribution function for the parameter a. When a member of the ith dynasty 
is born, he gives Xa a prior distribution^. We make the reasonable modelling 
assumption that this child's prior for Xa is Normal with mean oti and precision^!! e. 
The agent then updates his prior according to his observation of (X s ) t ^ <s<t , where 

t\ denotes the time of birth of the current child and t is the current time. 

If the agent knew the value of a, he would simply use a change of measure of the 
form (|2.2p . However, a is unknown, so the agent must weight each of the changes 
of measure according to his prior distribution for a. Hence at time t, agent i's law 
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This idea of dynasties has been used by iNakata ( 2007 ) , who considers an economy in which at 



any time point there are H young and H old agents. Each agent lives for 2 periods. Young agent 
h G {1, - -jH} has the same preferences and beliefs as the old agent h. He then considers a Rational 
Beliefs Equilibrium as explained by Kurz. However, all agents in his model live for exactly two units of 
time, in contrast to our ass u mptio ns. 

5 See lRogers fc Williams! (l2000h . IV. 38 for an account 



6 This is equivalent to having a prior distribution for a, since A is known. 
7 Equivalently, the prior has variance e _1 

8 Hence, all members of dynasty i begin life with the same prior 
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for the path has density with respect to the reference measure given by: 

Aj = ^exp ( - e -(\a - a*) 2 + \a(W t - W t t) - \{\af(t - t£))d(Aa) 

( (AW) 2 + 2c/eAVF - e(a^At \ 
ex P ( n? I a7\ ) 



e + At 2(e + At) 

where At and AT^ are given by: 

At = t-t{ AW = W t -W t i a = Xa (2.4) 



2.3 Deriving the State Price Density 

Associated with agent (or dynasty) i is a utility function, which we take to be 
CARA: Ui(t,x) = — ^-e~ liX e~ pt . Here, p is the discount factor, assumed to be the 
same for all agents. The agents seek to maximise the expected discounted utility 
of their consumption. Thus, agent i's objective is: 



max E 



to 



u t (t,4)K (2-5) 



where to is some start value, which we will later allow to go to — oo. K\ is the 
density derived in (|2.3|) . which jumps at each of the times T^. 

Th e objectives of the agents have the same form as the previous iBrown fc Rogers 
so its theory can be used to derive a state price density. In particular, by 
looking at the price of an arbitrary contingent claim we can deduce that: 

C s ^ = Ui(s, cj)A* 



where Vi is some Tt random variable! lj. Recalling our expression for Ui and 
taking logs, we obtain: 

togCt + !og^ _ _ j + log A f / 2 g x 

li li H * H 

Summing (|2.6p over i and using market clearing gives: 



lv^l l^log!/j loP* <5, 1 ^ log A 



2.4 A continuum of agents 

Recall that there are N different agents in our model. We will now let N tend to 
infinity so that we can examine the case in which there is a continuum of agents. 
We assume that ^ — has a finite limit and denote this limit by: 

r-i=ii m !yi 

N N ^— ' 7j 

9 We will shortly let to tend to negative infinity and when this occurs, the Tt Q will be trivial, thus i>i 
will just be a constant 

10 U' denotes the derivative of U with respect to its second argument 
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Abusing notation slightly, we use dj to denote the limjy_*oo 77- Hence: 

log C t + G' = -pt - T( ai X t + a 2 X 2 t ) + T Urn V J- log A} (2.7) 

N^oo ' — • Jv'yi 

where G" is some -measurable function. We now let to tend to negative infinity; 
Tt then becomes trivial, so G becomes a simple constanj^l. 

Only the last term in <\2.7h requires further development. Writing u 1 for the 
time since the the last person died in the ith dynasty, we obtain: 



T lim V — log Aj = T lim V 



1 rl l 



2 \e + u z 



( (W t -W t _ ui ) 2 + 2a*e(W t -W t _ ui )-e 
V 2(e + u*) 



(2i 



We assume that the mean of a 1 is given by (a) and further that the distribution 
of it 1 , a 1 and 7, are all independent. We further make the assumption that u has a 
density <£>(•), given by: 



9?(u) = A(e + it)Ae" 



where A = is chosen so that J °° tp(u)du = 1. Since </?(u) represents the 

probability of someone who is currently alive having age u, it follows that <£>(•) 
must be decreasing. This gives the inequality Ae > 1. 
Using our expression for ip, equation (12. 8p becomes: 

log (t = -G- Y(aiX t + a 2 X 2 ) -pt + \l ^ Wt ~ Wt ~^ ^{u)du 

2 J e + u 

+ (a)e / ^ L (p(u)du 

where G is some new constant. This then gives us: 

A 



where 



logO = ~G- T{aiX t + a 2 X 2 ) -pt + ^r, t + (a)eA& 



JO 

r?t = / (W t - W t _ u ) 2 Ae- Au d« 



By rearrangement and use of Fubini (see appendix), we are able to show that: 

Zt = x t 

Vt = X 2 t + e- xt [ \e Xs X 2 s ds 



11 We note that as to — > 00, the expression on the right of (|2.7|) is almost surely finite, so the left hand 
side must be as well. Since our Q and (fi)i<i<jv were only chosen up to a multiplicative constant, we 
may choose them to depend on to in such a way that as to —> 00 both £ and G' are a.s. finite. 
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Our final expression for the state price density is then given by: 

logO = —G — T( ai X t + a 2 X 2 ) -pt + ^ [TO 2 + e- Xt f \e Xs X 2 s ds] + (a)eAX t 

(2.9) 

= -G + BX t + CXf + U t -pt (2.10) 

where: 

B = (a)eA - Ten C = — - Ta 2 

and 

1 f* 
C/ t = -^e~ A ' / Xe Xs X 2 ds 

2 j— o 



-oo 



3 Asset Prices 

3.1 The Interest Rate Process 

We will use our state price density to derive the interest rate process. From Ito's 
formula, we have: 

^ = (B + 2CX t )dW t 

+ (C + ^-Xf -XU t -p- BXX t - 2\CX 2 t + \{B + 2CX t f)dt 

= [{-p + C+ l -B 2 ) + (-XB + 2CB)X t + (-2AC + ^ + 2C 2 )X t 2 - XU t ] dt 

where the symbol = signifies that the two sides differ by a local martingale. The 
interest rate is equal to minus the coefficient of dt in the above expansion, hence: 

r t = r(X t , U t ) = {p-C- X -B 2 ) + B(X - 2C)X t + (2AC - ^ - 2C 2 )X 2 t + XU t 

(3.1) 

Thus, our model gives us an interest rate process of the form: 

r t = a Q + aiX t + a 2 X 2 + XU t 

for some constants i = 0, 1, 2. Note that the interest rate process will depend on 
the behaviour of the dividend process in the past (via Ut) as well as on the current 
value of the dividend process. We therefore see that in some sense, high historical 
volatility generates high values of the riskless rate. 

3.2 The Stock Price 

We will now calculate the stock price. We have: 

s t =E° t \ r^du 



t Ct 

1 f°° r i 

- / E° ( U 5 U du (3.2) 

5* Jt L J 
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3.2.1 A PDE for the stock price 

Prom the form of Q and the Markovian structure, we will have that: 

( t S t = ( t h(X t ,U t ) (3.3) 

for some function h. This function will satisfy a PDE which we may determine by 
by observing that QSt + Jq Cs^sds is a martingale and applying Ito's formula. After 
a few calculations, we obtain the PDE: 

1 A 
= -h xx + (B + (2C - \)x)h x + Ki^x 2 - u )hu - r(x, u)h + (a + a\x + a 2 x 2 ) 

(3.4) 

Unfortunately, it does not appear to be possible to solve this equation in closed 
form, so we will resort to another approach. However, before we do this, let us look 
at some of the consequences of (|3.3p and (|3.4p . Suppose that under the real- world 
probability, P*, the OU process reverts to level a*, then we have that: 

dS t = h x dW t * + h x X(a* - X t )dt + h u (^-Xf - XU t )dt + ~h xx dt 



where W* denotes a Brownian motion under measure P*. After using (13. 4D we get 
that: 

dS t = h x dW t * + h x ({Xa* - B) - 2CX t )dt + r{X u U t )hdt - (a + a\X t + a 2 Xf)dt 
Hence, we see that the volatility and drift of the stock price are given by: 
_ h x (X t ,U t ) 

St " HXM) (3 - 5) 

- rjXt, Ut)h(X t , U t ) - (a + ai X t + a 2 A f 2 ) + (Xa* - 2CX t - B)h x (X t , U t ) 

h(X t ,U t ) 

(3.6) 



/'f 



We shall use these expressions later. 

3.2.2 Calculation of stock price via computation of conditional ex- 
pectation 

We will now proceed to determine the stock price via another method. Substituting 
the state price density from (|2.10p into (13. 2p . we obtain: 

S t = exp{-BX t - CX} -U t + pt} 

/•oo 

J E° [(oo + oi X T + a 2 X\) exp{BX T + CX\ + U T - pT}] dT 

On first sight it may appear that it is very difficult to get any further with this 
expression. However, if we can calculate: 

V T (t,X t ;6) :=E°[exp{0(a o + aiX T + a 2 x!) +BX T + CX% + J ^Ae A(s ~ T) X^ds}] 
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then we may differentiate with respect to to and set = to give: 

/OO Q 
exp{(e- A ( T -') - 1)17* - p(T - t)} — \ e=0 V T (t, X t ; 0)dT 



it 

We also define r = T — t. We will show that: 



V T (t, X t ; 0) = exp{la(r)X t 2 + b(r)X t + c(r)} 

where a, b and c are functions which we will shortly deduce. To deduce these 
functions, we will use a martingale argument. For t < T we define: 

M.J = E° [exp{#(a + aiX T + a 2 X|) + BX T + CX| + / ^Ae A(s " T) X 2 s ds}] 

= V T {t,X t -0)e^{f ^Ae^X 2 ^} 

J-oo * 

Now apply Ito's formula: 

/t 4 1 \ A 

-\e x(s - T ^X 2 s ds} [V t dt + V x dX t + -V xx dX t dX t + _ e A (*- r )^ 2 y^] 



= Mj xfdt - (\d{r)X? + 6(r)X t + c(r))cft 

+ (a(T)X t + 6(r))(dWi - AX t cft) + ^(a(r) + (a(r)JQ + 6(r)) 2 )dt] 

But (M t T )t<T is a martingale under Pq, so the coefficient of in the above expres- 
sion must be zero. Thus we obtain: 

L ^ -At x 1 2 

-a = — e Ar - Aa + -a 2 
2 2 2 

6 = a& — Xb 
c = \(a + b 2 ) 

The boundary conditions are given by: 

o(0) = 2(C + 0a 2 ) b(O)=B + ai c(0) = 0a o 

3.2.3 Solving the ODEs 

We now solve the ODEs. The first equation is a Riccati equation, so in order to 
solve we make the usual substitution: 

Substituting this into the ODE for a gives: 

1 ■• x • * A Xr 

^9 + ^9 + ~Y e 9 = 
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and the boundary condition becomes: 

-5(0) = 2(C + 9a 2 )g(0) 
We can solve this equation using Maple to obtain: 

g (u) = e~ Xu \(VXAY 1 {2^/aJX) - 2(C + 9a 2 )Y 2 (2^/A/X)) J 2 {2e~ Xu l 2 ^/A/X) 



- (VXAJ^y/A/X) - 2{C + 9a 2 )J 2 (2^/A/X))Y 2 (2e' Xu / 2 ^/A/X) 

where J{ and Yi are Bessel functions of order i of the first and second kind respec- 
tively. Turning now to the ODE for 6, we may use our solution for a to deduce: 

b+ h + Xb = 



Rearranging gives: 



which we can solve subject to 6(0) = B + 9a\ to give: 

e XT g{r) 

Finally, we obtain: 



c(r) = 9a + j^(a(T') + b(T') 2 )dT' 



Thus we have completely solved the ODEs. In order to calculate the stock price, 
we need to find We therefore need: 

d A = e -A«[_ 2a 2 Y 2 (2^/A/X)J 2 (2e- Xu / 2 ^/A/X) + 2a 2 J 2 (2^/A/X)Y 2 (2e' Xu / 2 y^A/X)] 
and also: 

|| = -A^+ e - Au [-2a 2 y 2 (2yi7A)(AJ 2 (2yi7A e - Au / 2 )-v / AAe- Au / 2 J 1 (27I7Ae- Au / 2 )) 
+ 2a 2 J 2 {2^A/\){XY 2 (2^A/\e- Xu ' 2 ) - VAXe-^Y^^/A/Xe-^ 2 ))] 

We may then calculate expressions for First note that: 

dV ,19a n 9b 9c x A . . o , / ^ / m 
= (o + 7ti X * + 7S») exp{-a(r)X 2 + 6(r)X t + c(r)} 



But: 



90 l 29« 1 9« ' d9> ^ L 2 

9c . , f T 1 ,da . ,, ,.96, M , , 

09^ = ao + J 2^ (r) + 26(T) ^ (r) ) dT 

9& M= g(0) (fl + 0ai)gg(O) (g + gai) g(0) 9g 
90 lTj ai e A - 5 (r) e A - 5 (r) e A - <?(r) 2 90^ 

da (r) = ^ (T) I g(r) —(t) 
d9 { ' g(r) ^ g(ryd9 { ' 
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So finally we have: 



poo 

S t = exp{-BX t - CX?} \ exp{-pr - (1 

Jo 



-At 



Wt} 



^\%^ + te Xt + ex p4 a(r)x * + b{T)Xt + c(T)}dT (3 - 7) 



30 ' 



This is as far as we can get with the expression for the stock price. We see that 
the stock price depends not only on the dividend at time t, but also on Ut, a term 
reflecting the behaviour of (X s )_ oc < s <t. This is as we would expect, since agents 
need to use information from the whole of their lifetimes to make better estimates 
of the mean to which X is reverting. From properties of the OU process, we see 
that if X t reverts to mean a then, since X is stationary, we have X t ~ N(a, ^). 
Hence, 

This indicates a sensible value for Ut, which will be helpful for when we begin to 
look at numerical examples later on. 

3.3 The Bond Price 

The time-t price of a zero-coupon bond which has unit payoff at time T is given by: 



Ct \T 



exp 



BX t - CX 2 t -U t (l- e 



-\t\ 



pr 



V T (t,X t ;6 = 0) 



Using our expression for V T (t,X t ;6 = 0), we obtain: 



exp 



(-a(r) - C)Xf + (6(r) - B)X t + c(r) - pr - (1 - e~^)U t (3.8) 



where the functions a, b and c are all evaluated using 8 = 0. 



3.4 Remarks on the case in which a is known 

Note that if we let e — > oo, then this corresponds to the case in which all the agents 
are certain that they know the value of a. By taking the limit in our expressions for 
the stock price, bond price and riskless rate, we can deduce expressions for these 
quantities in this limit. We note further that if the agents are sure about the value 
of a and this value corresponds to the true value, a* , then the expressions we obtain 
will be the same as those for the model in which the true value of a was known to 
all the agents. 



4 Numerical Results 
4.1 Calibrating the model 

We now investigate the qualitative behaviour of the stock price as we vary the 
parameters in our model. To do this it is important that we choose a suitable 
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region over which the parameters can vary. We will restrict to the case in which 
ao = a\ = 0, so that we have simply St = a^X 2 . This ensures that the dividend 
process remains positive. Note further that the state price density (|2.10|) only 
depends on the product Ta>2 rather than the individual T and 02- Although the 
dividend process does depend on 02, changing 02 simply corresponds to the changing 
the units in which we measure the dividend process. Hence, we may choose 02 = 1. 

Some of the parameters are relatively easy to choose, such as A and p, for which 
we choose A = 2 and p = 0.04. However, other parameters, such as T are much 
harder to determine. We are only interested in ensuring that the parameters are 
of the correct order. In order to do this, we will begin by assuming that (a) = a. 
Furthermore, for the purposes of calibration, we will consider the limit as e —* 00, 
which corresponds to the case in which agents are sure that they know the true 
value of a. This leaves the parameters a and T which we still need to determine. 

One way to determine these parameters would be to choose them in order to 
match various moments from empiri cal data, such as t h e mea n price-dividend ratio; 
this was the strategy employed in iBrown fc Rogers! (|2009h when we considered 



the equity premium puzzle. Ideally, we would use the same method here, but 
unfortunately our stock price is much more complicated. Thus, computing a given 
stock price requires the numerical computation of an integral. To work out the 
mean price dividend ratio, we would then need to compute a further integral as 
we averaged over the values of the driving Brownian motion. We would then vary 
the parameters and calculate the expected price dividend ratio each time in an 
attempt to find a realistic set of parameters. Given the additional complexity of 
this problem and the fact that we are only interested in determining parameters 
that are of the correct order, we will proceed in a different manner. 

We first note that the interest rate process has a particularly simple form, which 
we can use to get a simple expression for the expected riskless rate. We can match 
this with the mean riskless rate from the Shiller data set. 

Note that we are considering the case in which ao = a\ = 0, a2 = 1 and the limit 
as e — > 00 and hence Ae — ► 1, B — ► (a) = a, C — > —T. Substituting into expression 



(|3.ip gives: 

r = (p + r - \a 2 ) + a(A + 2T)X t - 2T{\ + T)X 2 

Thus, the expected riskless rate is given by: 

Er = (p + T - ±a 2 ) + a 2 (A + 2r) - 2r(A + T)(a 2 + 

To determine T, we compare a CRRA agent (where we know a reasonable value 
for the constant of relative risk aversion) with a CARA agent. If we consider a 
single agent model in which the value of a is known, the stock price will be given 
by: 

r°°u'(s t ) 
So = E J uW) 5tdt 

Since we just want our parameters to be of the correct order, it is sufficient to check 
that the behaviour of 

wW) 6t (41) 
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in the neighbourhood of X = a is the same for both the CRRA and CARA case. 
If we set Xq = Xf = a then clearly (|4.ip will be the same in both the CRRA and 
CARA case. We therefore impose the requirement that a small change in Xt from 
Xt = a has the same effect in both cases, leading to the condition: 

U'cRRAi " 2 ) _ U'cARA( a2 ) 
UcRRA( a2 ) U' CARA {a 2 ) 

Which leads us to the condition: 

r = 4 ( 4 - 2 ) 

Since we know a sensible value for R is R = 2, this equation gives us an equation 
from which we can determine T and a. Substituting in our expression for the riskless 
rate yields the cubic equation: 

r 3 

K r ) = -r + 2Rv2 + ( Er ~ P + 2R ( X ~ !)) r + \ R ~ RX = 

A 

We will choose R = 2. We also choose Er = 0.01, as given by the Shiller data set. 
We may then note that 1(0) < and jj^ > for T > 0, hence there is a unique 
positive solution to the above equation, which we can easily compute. Computation 
shows that the correct V to choose is T = 0.49 which we take as our default value. 
This gives a=2.01. 

To summarise, the default parameters we choose are: ao = ai = 0; 02 = 1; A = 
2-p = 0.04; e = 1.0; T = 0.49; (a) = a = 2.01. We also choose X t = a, U t = 

— (a 2 + ^). We then vary the parameters and examine the behaviour. 

4.2 Comments on Results 

Figure [1] shows that the stock price is decreasing in A. Recall that A is the parameter 
which tells us how quickly the dividend process returns to its mean. Hence, a lower 
value of A means that the dividend process is more likely to reach high values, so 
is worth more to the agents. However, A is also a parameter used in specifying 
the distribution of the lifetime of the agents. Increasing A therefore decreases the 
expected lifetime of the agents. Each child in the dynasty therefore has less time 
to learn about the unknown parameter a and this increased uncertainty amongst 
the agents also means that the stock price decreases as A increases. 

Figure [2] shows that as e increases, so does the stock price, which is to be 
expected since if the agents know more about the dividend process (i.e. their 
beliefs have a higher precision), the stock should be worth more to them. Similarly, 
Figure [3] shows that the larger the value of p, the less the stock is worth. A large 
p indicates that the agents are impatient and want to consume their wealth in the 
near future, making the stock less attractive. 

Figure 0] exhibits the dependence of the stock price on (a). Recall that Xt and 
Ut are kept fixed as we vary (a). A small (a) indicates that the agents think the 
level to which X reverts is low. Thus, since we do not change Xt, a low value of (a) 
relative to X indicates that X is currently abnormally high and so the dividends 
are abnormally high. Thus, the agents are keen to hold this stock. Furthermore, 
the relatively high level of X means that the agents have a large amount of dividend 
with which to buy the stock. 



14 



Figure [5] may at first seem surprising, since it shows that the stock price is 
increasing in the risk aversion, T. However, we recall that all agents have a CARA 
utility and furthermore, the parameters of our model are chosen so that the dividend 
process is non-negative. On the one hand, a larger value of T means that the value 
of the dividend process becoming larger are valued more highly than before. The 
downside of holding the stock is limited, since the dividend process is always non- 
negative. This explains the behaviour shown in[5J 

The volatility surface^! in Figure [6] shows that the volatility appears to be 
increasing in both Xt and Ut- This seems reasonable: if the dividend process has 
been varying greatly in the past, then Ut will be large, and in this case we would 
expect the stock to have a larger volatility. Also, if Xt is small and Ut is large, 
it means that either the dividend has been varying greatly, or the value of Xt 
is abnormally small, so in these large volatility should not be surprising. 

However, increasing Xt means that Ut no longer implies that the dividend process 
has been varying so much - it just tells us that Xt is typically large. This explains 
why increasing Xt will decrease the volatility. 

5 Conclusions 

We have introduced a new model in which the dividend of the stock obeys an OU 
process for which none of the agents know the mean. We derived a state price 
density and were able to use this to price the stock and a bond. We also were 
able to deduce an interest rate model. We produced graphs which illustrated the 
dependence of the stock price on the various parameters. The behaviour shown 
in these graphs seemed very reasonable. We also looked at how the parameter 
certainty case could be viewed as a special limit of the parameter uncertainty case. 

Extensions to this work include using a different utility function for the agents; 
a CRRA utility would be a natural choice. In section 12.41 we also had to assume a 
quite specific form for the distribution of the lifetimes of the agents. An obvious 
improvement would be to consider the problem with a different distribution of 
lifetimes, in particular one that did not depend on the parameters of the dividend 
process. Unfortunately both these generalisations appear to make the calculations 
intractable. 
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APPENDIX 

A Stochastic Integrals 
A.l Calculating £ 4 

Recall that £t is given by: 

pOO 

6 = / (W t - W t . u )\e- Xu du 
Jo 

By change of variables, 



tt = W t - eT xt / \e Xs W s ds 



So substituting from (12. ip gives: 



& = W t + X - e~ xt 



t rt rs 

X s Xe Xs ds + / Ae As / \X r drds 

-oo J-oo jO 



(A.l) 



But the final term in the above expression is: 

f xeXs F 

J -oo JO 



-\t f Xe Xs f xx r drds 



r"- v / / Xe Xs XX r drds - e~ xt [ [ \e Xs XX r drds 



s=—ooJr=s J s=0 J r=0 



Applying Fubini, we obtain: 

pO pr pt pt 

e~ xt / / Xe Xs XX r dsdr - e~ xt / / 

Jr=—ooJs=—oo J r=0 J s=r 

Computing the integral with respect to s gives: 



e -xt 



[0 ft ft 

[ / Xe Xr X r dr - e xt / X r Xdr + / Xe Xr X r dr] 

J -oo Jo Jo 



e xt i Xe Xr X r dr- [ XX r dr 



Substituting this into (|A.1|) gives: 

= W t + X - [ XX r dr 

But recalling ()2.ip . we obtain: 

tt = X t 



o 
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A. 2 Calculating rjt 

Recall that rjt is given by: 



poo 

m= / (W t -W t - U f\e- Xu du 
Jo 



Changing variables we obtain: 



Vt = e~ Xt 



f (W t - W r ) 2 \e Xr dr 

J — oo 



Substituting from (|2.ip gives: 

[{X t -X r )+ I XX s ds] 2 Xe Xr dr 

-oo J r 



Vt = e~ Xi 



e- xt 



f {X t - X r f\e Xr dr + 2e~ xt f (X t - X r ){ f XX s ds)Xe Xr dr 

/t rt 
( / \X s dsf\e Xr dr (A 
-oo J r 



We will now apply Fubini to two of these terms to deduce an expression for 
Firstly, we work on: 



r=— oo J s=r 



X t I XX K ds\e Xr dr 



By applying Fubini, we obtain: 



f X t X s f X 2 e Xr drds 

J s=— oo J r=— oo 



' s=— oo j r=— oo 

ft 

X t X s \e Xs ds 



Putting this into (|A.2h gives: 



m = X 2 t + e~ xt / Xe Xr X 2 dr 



2e~ xt f X r ( f XX s ds)Xe Xr dr 

J —oo J r 

+ e' xt I j XX s ds I XX v dvXe Xr dr (A 

J —oo J r J r 



The final term is: 



ft rt rt 

2e~ xt I / / XX,XX v Xe Xr dvdsdr 



r=— oo J s=r Jv=s 
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where we have halved the area of integration in the dvds integral. Applying Fubini 
yields: 



2e 



-At 



\X,\X v \e Xr drdvds 



s=— oo J v=s Jr=—oo 
t ft 

2r _A ' / / \X s XX v e Xs dvds 

=— oo J v=s 
t r-t 

2<r xt I XX r e Xr / \X s dsdr 

T— — oo J s=r 



Substituting this into ()A.3|) gives: 



Vt = X 2 t + e- xt f \e Xs X 2 s ds 

J — oo 
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Figure 1: Graph of S t against A. 
Stock Price vs lambda 
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Figure 2: Graph of S t against e. 
Stock Price vs epsilon 
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Figure 3: Graph of S t against p. 
Stock Price vs rho 
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Figure 4: Graph of S t against (a). 
Stock Price vs alpha 
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Figure 5: Graph of S t against V. 
Stock Price vs Gamma 




Gamma 
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Figure 6: Volatility Surface 
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